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Abstract 



We discuss two technical issues related to the gap equation in high-density 
QCD: i) how to obtain the asymptotic solution with well controlled approx- 
imations, and ii) the renormalization of four-quark operators in the high- 
density effective field theory. 



Recently Son obtained the leading exponential behavior of the superconducting gap in 
QCD at asymptotically high density using both an indirect renormalization group argument 
and a direct QCD calculation Subsequently many papers confirmed Son's result [@-|J. 
However, to the best of our knowledge, the analytical determinations of the gap suffer from 
two flaws. First, the gap equation is divergent and so must be regularized and renormal- 
ized. Most treatments have taken the baryon density as a sharp cutoff. This might cause 
some concern since it leads to the possibility of contaminating the low energy physics of 
the gap with ad hoc high energy physics. Issues of cutoff sensitivity have been addressed in 
Ref. || and Ref. ||. Second, the solution of the gap equation at momenta large compared 
to the gap has been obtained by assuming (what appears to be) a particularly unhealthy 
approximation which allows the integral equation for the gap to be expressed as a simple 
differential equation. The goals of this paper are modest. Using cutoff regularization we 
define a renormalized gap equation and we find the exact asymptotic solution, thus excising 
the flaws contained in previous determinations. Our results confirm Son's original analy- 
sis. We also obtain the asymptotic solution using dimensional regularization with minimal 
subtraction. 

Many degrees of freedom, including antiparticles and hard gluons, are not dynamical on 
the Fermi surface. Hence it is sensible to work with an effective theory of QCD appropriate 
to the scales in question. Explicit construction of the effective field theory appropriate 
for momentum scales below 2/x, can be found in papers by Hong ||. The fermions in this 
effective theory live on the two-dimensional Fermi surface and so depend only on the parallel 
momentum, q». The gluons on the other hand propagate in directions perpendicular to the 
Fermi surface as well and therefore also depend on the perpedicular momentum, q±. Hence 
we should treat the effective field theory as a superposition of two-dimensional theories, one 
for each direction on the Fermi surface, interacting through four- dimensional gluons and 
contact operators. Only the graphs shown in Fig. (|l]a) and (|l]b) need be calculated if we are 
interested in the leading exponential behavior of the gap. The sum of these graphs is (after 
rotating to Euclidean space) 

r d\ A(g„) Wf£q±( 1 , _J__] , fX m 

[PV J (27r)^f+A(g||)H 3 J (2ir) 2 \ql + f Mj|p - 9o|/|<zi \ Q± + M d) J' 1 J 

where the propagators within the parentheses represent magnetic and electric gluon ex- 
changes, respectively, and 

Ml = ^ (2) 

The effects of hard gluon exchange, antiparticle exchange and any residual gauge dependence 
are represented by the coefficient, D, of a four-fermi operator in the two-dimensional effective 
theory. 

It is clear that the integration over q± is divergent. We will first regulate the gap equation 
with a sharp cutoff, Aj_. In principle there is a cutoff associated with the parallel integration 
as well and therefore in general we have D = Z?(A_l, An). It is straightforward to do the 
integration over q± and q\\. We obtain 

Afo.) = / A " dg *«>] ( log + W A,,)) (3) 
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FIG. 1. T/ie leading diagrams contributing to the gap equation. The solid square denotes the 

gap, A, while the solid circle denotes an insertion of the counterterm, D. 
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Since A(po) is independent of the cutoff, D runs according to 

L»(A ± ,A||) = J D(? 7) A||)-61og^. (5) 
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Naive dimensional analysis suggests that D(A± = 2/x, An) is small and can therefore be 
dropped at leading order. The effect of the running of D appears at next order ||]. We 
then have 
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D(2/i,A||) = D(A|,). 



2/i iv//y ' 

Dropping all subscripts for simplicity the gap equation becomes 



A(p) = ^C 2 / A rfg 



A(g) 



g 2 + A(g) 2 V \P 2 -Q 



M 2 \ 
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(7) 



Because of the log singularity at g ~ p, the integral is dominated by momenta q ~ p. 
Therefore in the asymptotic region defined by p ^> A(p), we can take g ^> A(g) under the 
integral. Hence asymptotically the gap equation becomes the homogeneous integral equation 

A c)4 c2 r7 Aw ( ios p^ +c(A) )' <8) 

Note that although the integral equation is homogeneous, we have by necessity introduced 
an infrared cutoff, A, which we will see is related to the overall normalization of the gap. In 
order to proceed^, consider the derivative of the gap equation: 



1 The usual way of proceeding is to make the approximation log \p 2 — q 2 \ 
log q 2 9(q 2 — p 2 ). We do NOT make this approximation in this paper. 



logp 2 6>(p 2 — q 2 ) + 
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A'(p)=C 2 pf 



A(g) 
A q (q 2 — p 2 



(9) 



Note that since the counterterm is no longer present, this integral is convergent as A — > oo. 

We make an ansatz of the form A(p) = p z with z a complex number. Inserting this 
solution in eq. (H) leads to the indicial equation: 



z = C 2 p 2 -* dqj^-^. 

J A {q z — p z ) 

The integral is straightforward to evaluate. It is given by 
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where the contour in the complex q-plane is taken to enclose the poles at p — ie and — p + ie 
while avoiding the branch point at the origin. For A <C p C A we can take the limits A — ► 
and A — > oo since the corrections are suppressed by powers of p/A and A/ p. The resulting 
indicial equation is transcendental 



-C 2 cot—. 



(12) 



Therefore the exact asymptotic solution of the gap equation is p z with z satisfying eq. fH2| 
For small z this has the solution 



2! = ±iC. 



(13) 



The general asymptotic solution to the gap equation in the weak coupling limit can then be 
written as 



A* 

A(p) = A sin (C log — ), 
p 



(14) 



where the constants A and A* are to be determined. 

We will now determine the constant A. In the region A < p < 2/i, A(p) has one 
maximum, which we will assume is at p = A with A < A < A. This determines A = A(A). 
We can find A by plugging the general solution, eq. ([14]) , into eq. @ which leads to the 
equation 



= A'(p)L A = C 2 AA(A) 



A straightforward computation then gives 
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where 
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A* 

x A = Clog— ; 



A* 

X A = Clog— . 



(17) 



The expression under the sum is exponentially suppressed in the QCD coupling g. Therefore, 
to leading order xa = 7r/2 and we can identify A with A. This determines A = A(A) and 
the asymptotic solution is 



A* 

A(p) = A(A)sin (Clog—); 

V 



A(A) = A = A*exp(--). 
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We now determine A*. We can rewrite eq. (|S|) as 



where 



(18) 



(19) 



X 2 = .M 2 exp( J D(.M)). 



(20) 



The condition that A(p) be independent of the choice of cutoff leads to the renormalization 
group equation 



A-^- (D(A) - D(M)) = -C tan(C log ^~ 



dA 
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D(A)-D(M))+log 



A 2 



(21) 



where we have used the asymptotic solution, eq. fll8l). The solution of this equation is 



D(A)-D(M) = - 



tan(Clog^- Sm(C1 ° g ^' 



A^ cos(Clogf) 



- log 



M 2 
A 2 



(22) 



Notice that the difference between counterterms evaluated at different choices of scale is 
of order g 2 . The counterterm has a strong cutoff dependence except when A ~ M.. The 
peculiar running of the counterterm is very similar to the running of a counterterm which 
arises at leading order in effective field theory treatments of the three-body problem in 
nuclear physics 0. 

Choosing A = Ai, the renormalized gap equation is 



A(p) 



-c- 



M 



dq 



A(q) 



g 2 + A( q y 



log 



p 2 — q 2 \ 



(23) 



We can now find A* by plugging the general solution, eq. (|T8|), into the asymptotic form of 
the renormalized gap equation, eq. (|23|) , evaluated at p = A*. We then have 



1 rMrja A* 

= A(A*) = -C 2 A(A) / ^ sin (Clog—) log 

2 Ja q q 



|A* 2 -g 2 | 



(24) 



A straightforward computation then gives 
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where 
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The expression under the sum is exponentially suppressed in the QCD coupling g. Therefore, 
to leading order Xj^ = and we can identify A* = M.. 
The final form for the asymptotic solution is thus 



M 

A(p) = A(A)sin (C log — ) 

V 



and the gap is 
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(27) 
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This is Son's result aside from the additional contribution to the prefactor from the 
counterterm, which is expected to be a number of order one. 

One may wonder whether use of cutoff regularization in dense QCD violates important 
symmetries. Consistent implementation of dimensional regularization would erase these 
concerns. We will see that dimensional regularization with minimal subtraction is a quick 
way of obtaining the asymptotic solution directly from the gap equation. We can continue 
the four-dimensional measure to an 2 + n-dimensional measure 



d 2 qn f d n q±_ 



(2vr) 2 J (2?r) n ' 

The gap equation relevant to dimensional regularization with minimal subtraction is 



(29) 



A(P||) 
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(27r) 2 gf + A(g|| 
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The integrals of the gluon propagators in n-dimensions are 



d n q ± 1 1 
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(31) 



where e = (2 — n)/3 and A is a renormalization scale. Absorbing the 1/e pole into the 
counterterm, we can then define the regularized gap equation 



A(p) = \C 2 



dq 
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(32) 
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where 
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This equation was obtained in Ref. [|]]. The counterterm runs according to 



D MS (X) = D MS ( V )-Q\og^. 

T] 2 



(33) 



(34) 



Physical quantities are A-independent so we choose A = 2/i. We again consider the asymp 
totic gap equation 

1 



(35) 



q \ \p — q 

Say the asymptotic solution is of the form A(p) = p z . All scales then appear multiplied 
by power law divergences which vanish in minimal subtraction. Hence it is appropriate to 



return to the unsubtracted expression for the log in eq. (31) . We then obtain 

C 2 T(e] 



V 



Q 



d 1-( 2 
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Again this integral is straightforward to evaluate and leads to 

r ( e -f) r (f) 



p z = c 2 r( e 

which as expected reduces to 



p z 26 . z, 

cos7r(e , 
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in the limit e — ► 0. Notice that the original T(e) pole from the integration over perpendicular 
momenta cancels a 1/T(e) zero from the integration over q. The asymptotic solution is again 



A* 

A(p) = A(A) sin (C log — ) 

P 



(39) 



where the prefactor is fixed by the argument presented above. Now fixing A* is trivial since 
there is only one scale in the problem. We have 



A* = A^exp 



D MS (2fi) 



from which follows 



A 
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(40) 
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Evidently the counterterm does not depend on the renormalization scale in the longitudinal 
direction in minimal subtraction. As argued previously, the counterterm is of order g° and 
can be dropped at this order in the expansion. 
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